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By considering mirror oscillation in a "tripod-scheme" laser-atom system, we advocate explorative 
studies of driven Dirac-like equations. Both analytical and numerical studies show that mirror oscil- 
lation can be used to drive an effective spin-orbit interaction and hence control the amplitude, the 
frequency, and the damping of the cold-atom Zitterbewegung oscillation. Our results demonstrate 
an interesting coupling between the mirror mechanical motion and a fundamental quantum coherent 
oscillation, opening up new means of matter wave manipulation. 

PACS numbers: 03.75.-b, 32.80.Qk, 71.70.Ej, 37.10.Vz 



I. INTRODUCTION 

4^* ■ 

^ ' The jittering motion of a free relativistic electron predicted by the Dirac equation, called Zitterbewegung oscillation 
^ ■ (ZB) Jj, is a truly fundamental quantum coherence effect. Directly observing the free-electron ZB is, however, 
^H' practically impossible due to its extremely high frequency and small amplitude. For this very reason, studies of 
I— I, ZB effect necessarily call for quantum simulations based on Dirac-like equations with (effective) spin-orbit coupling, 
including those involving band electrons in graphene |:2] , cavity electrodynamics [3| , single trapped ion 0, Q j as well 
\ as ultracold atoms 

• To actively explore ZB-related physics, it is necessary to go beyond passive simulations of the known ZB effect. 
— . ] Here we advocate to consider Dirac-like equations driven by an external field. This is feasible due to the precise 
controllability of laser-atom interaction. Given the current vast interest in dressed matter waves Q , this topic is also 
timely because it can reveal how a driven effective spin-orbit interaction may be used as a new means of matter-wave 
^f-^ manipulation. Furthermore, as proposed below, a driven Dirac-like equation can be achieved via mirror oscillation, 
' thus directly coupling the fascinating cold-atom ZB with the mirror mechanical motion. This interesting quantum- 
, classical interface may lead to a novel setup of optomechanical systems jlQ] . with spin-orbit interaction also included. 
In particular, to realize a driven Dirac-like equation we propose to add oscillating mirrors [see Fig. 1(b)] to a recent 
cold-atom-ZB scheme [3] that involves tripod-scheme cold atoms (TT| - [l3| interacting with three laser fields. This 
• ■ modification induces a coupling between the mirror mechanical motion and the cold-atom matter wave, yielding a 
' time-dependent effective spin-orbit interaction. It is shown that the amplitude of the cold-atom ZB can then be either 
^ . enhanced or weakened. Such control over the ZB amplitude is also related to a dynamical realization of an effective 
"spin-helix" Hamiltonian [l^ and the celebrated phenomenon of "coherent destruction of tunneling" (CDT) [l^ in 
driven systems. More remarkably, the quick damping of the cold-atom ZB, which hinders experimental studies, can 
also be dramatically suppressed. Extending the lifetime of cold-atom ZB might be useful for finding its applications in 
precision measurements or sensing. It should be also stressed that although we present our findings in the cold-atom 
context throughout this paper, they should have direct analogs in other alternative ZB realizations mentioned above. 



II. DRIVEN DIRAC-LIKE EQUATION 

Consider then the interaction of a tripod-scheme cold-atom interacting with three laser fields [ill, 113 [^^e Fig.l 
(a)]. The four internal atomic levels are denoted \n), with n = — 3, with the three states |1), |2), and |3) being 
degenerate magnetic sub-levels on the ground state. Each of the three laser fields has an appropriate polarization 
and induces a transition |0) ^ \n), with the Rabi frequency r2„, ri = 1 — 3. In the interaction picture and under 
the rotating wave approximation the internal Hamiltonian is given by -ffRWA,4 = ^'^n=ii^n\0) {n\ -f 17* |n)(0|). This 
Hamiltonian possesses two dark states \Di), I-D2), with zero eigenvalue and zero overlap with the excited state |0). 
For sufficiently large Rabi frequencies r2„ and for slow translation motion (for the system parameters associated with 
all our computational examples below, it can estimated that the characteristic Rabi frequency should be around 10^ 



2 




FIG. 1: (Color Online) (a) Level structure of tripod-scheme cold atoms, (b) Schematic setup involving oscillating mirrors and 
three laser fields interacting with cold atoms. As shown in the text, this leads to a driven Dirac-like equation. 



Hz or higher), the internal state can remain in a two-dimensional dark subspace spanned by \Di) and \D2). In the 
dark state representation which can be space-dependent, the translation motion of the atom effectively experiences a 
non-Abelian gauge field. The total stationary Hamiltonian then becomes 

Hs^:^{p-Af + ^ + V, (f) 

where m is atom mass, p is atom momentum. Ay = ih{Di\V Dj) is the effective non-Abelian 2x2 potential matrix, 
$ is the Born-Huang scalar potential, and V is the external potential. For convenience we assume ^ + V ~ constant 
and adopt the laser configuration studied in Ref. [l3|, with two laser fields counter-propagating along x, and another 
propagating along z, ili = 17o sin(^)e"*'='^/%/2, ^2 = rJo sin(^)e*'='^/V2, = f^o cos(Oe*'''^ and cos(C) = ^2-1, 
where ki is the laser wavevector. It then suffices to consider the x — z plane in Fig. 1(b) with two unit vectors Cx and 
Cz. As such, p = PxGx + Pz&z, and the Hamiltonian Hs becomes 



^S = 7i 1 [PxCTx+PzCrz), (2) 

2m TO 

where k, = {V2 — l)ki, and ax,z are the standard Pauli matrices. The second term in Eq. ([2]) represents an effective 
spin-orbit coupling that resembles the Dresselhaus spin-orbit Hamiltonian (other types of spin-orbit interaction, such 
as a mixture of Dresselhaus and Rashba coupling, can be obtained by considering different laser configurations [3]). 
When the effects of the p^/2to term are small as compared with the spin-orbit coupling term, the Schrodinger equation 
for Hs will be dominated by the second spin-orbit term in Eq. ^ , which is linear in px or pz . In this sense a Dirac-like 
equation (with two-component spinors) is obtained. 

To arrive at a driven Dirac-like equation, we consider two oscillating mirrors joined together, schematically shown 
in Fig. 1(b). If the mirrors are moving at a velocity v/2 = {v d/ 2) cos{u! dt) in the z direction, then effectively the fia 
field is moving towards the atom at a velocity v. As indicated below, such a movement is very slow as compared with 
the speed of light and hence relativistic Doppler effects can be safely neglected. In a moving frame where the laser 
source is stationary, Eq. ([2]) applies. We then make a Galileo transformation back to the laboratory frame, obtaining 
the following driven Hamiltonian 

p2 

-ffdrivon ^ t; ^ {[Pz " "^^'d COs(wrft)]crz -|- PxCTx}- (3) 

Zm m 

This driven Hamiltonian can be alternatively derived by considering non-Abelian geometric phases induced by the 
moving mirrors [T6| . When the spin-orbit interaction in Eq. ([Sj (which is now time-dependent) accounts for the main 
physics, the Schrodinger equation for i/drivcn simulates a driven Dirac-like equation. In the following two driven cases 
will be elaborated as examples. 
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FIG. 2: (Color Online) Effects of a high-frequency mirror oscillation on cold-atom ZB. The initial momentum po is taken in 
the z direction in (a) and in the x direction in (b). Insets are the typical wavepacket profile at early and later times, with the 
initial variance in x and z being around 10.0. Throughout t is in units of m/{hK^), x and z axe in units of u)d ~ 50 in 
units of VtK^ /m. 



III. CONTROLLED COLD-ATOM ZITTERBEWEGUNG 



Figure 2 depicts the time-dependence of the expectation value {x) or (z) of the cold atom, calculated from the 
evolution associated with -ffdrivcn- The initial state is given by |^'(0)) = (l/-\/2, i/v^)"^ / 'ip|p)G(p), where (r|p) ~ 
e'P '", and G(p) is a Gaussian distribution. This initial state represents a two-component spinor in the dark state 
representation times a Gaussian wavepacket. The shown oscillations in {x) or (z) without mirror oscillation (solid 
lines) are the expected cold-atom ZB effect. But remarkably, the ZB amplitude and frequency are seen to be strongly 
affected by introducing the mirror oscillation. In particular, in the case of Fig. 2(a), the initial wavepacket is set 
to move along the z direction, ZB occurs in (x), with its amplitude tuned down extensively and with its frequency 
unchanged. In the case of Fig. 2(b), the initial wavepacket is set to move along the x direction, ZB occurs in (z), with 
its amplitude tuned up extensively and at the same time its frequency much decreased. Note that in all the shown 
examples the ZB frequency is much smaller than the mirror oscillation frequency Wd- 

To understand the results in Fig. 2, we adopt a perspective afforded by a standard high-frequency approximation. 
That is, when 1 /w^ is much smaller than all other time scales, the driven dynamics can be approximately described by a 
static effective Hamiltonian ^f^fiven ^'^^ averaging of the time-dependent term of -ffdrivcn, i-e-, —hnvd cos{uJdt)az = 
i^drivcn - Hs = A{t). To the first order of l/w^, we obtain iT^^^^n = (wd/27r) J^"'^'^'' dt e'^^^*'' Hse-'^^*\ where 
F{t) = J^dt'A{t')/h [13. Specifically, 

2m m ujd 

where Jq is the ordinary Bessel function of order zero. Let 0{p) be the angle between the vector p = PxJo{^^^)^-x + 

PzCz and the x axis, and let /3(p) — 7r/4 — 6'(p)/2. The eigenstates of i?drivcn ^'"'^ then found to be 1?/;^) = 
[cos(/3), sin(/3)]'^|p) and IV'") = [sin(/3), — cos(/3)]'^|p), with their eigenvalues given by iJ±(p, p) = p^/2rn± ?iK|p|/m. 
Because each p component in the initial state |^(0)) can be expanded by the eigenstates IV"^), it can be predicted 
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from iJ^fiven t^^t total state at time t is given by 

ig+(p,p)t iE_(p,p)t _ "I 

X e « |^+) _ ie n IV' ) • (5) 

Fm^ther using r ~ ifiVp, one easily finds 

(r(t)) = *?i(M'(t)|Vp|*(t)) =r'' + p°Vm 

+ dp|G(p)HVp0(p)]{l-cos[..(p)i]}, (6) 

where the last oscillating term represents cold-atom ZB, with the angular frequency w(p) = [i?+(p, p) — i?_(p, p)]/?i = 
2k|p|/to, r° is the initial expectation value of r, and p° is the initial momentum of the Gaussian wavepacket. 
We now come back to the results in Fig. 2. In the case of Fig. 2(a), p^; sa 0, p « Pz^-z, then a;(p) sa i^(p) and 

rr / 2KVd \ 

Vp0(p) « ^^e,. (7) 

Pz 

Substituting these two relations into Eq. (jS]), one directly obtains that ZB is along x, its amplitude is proportional 
to the factor >Jo{^^^), and its frequency is independent of ojd- Quantitatively, for the two shown examples with 
2KVd/uJd = 1.52 or 2.22 in Fig. 2(a), the ZB amplitude should decrease by a factor of 2.0 or 10.0 as compared with 
that without mirror oscillation (v^ = 0), in excellent agreement with the numerics. Likewise, in the case of Fig. 2(b), 
0, p w pxix, leading to w(p) sa Jo(^^)w(p) and 
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PxJoi 



Vp^(P) « ^ w2..., e^- (8) 



Equation © then predicts that ZB is now along z, its amplitude is proportional to the factor l/Jo(^^^), and the 



ZB frequency is proportional to Joi^^^)- For the two examples with 2KVd/^d = 1-14 or 1.52 in Fig. 2(b), the 
ZB amplitude should be enhanced by a factor of 1.41 or 2.0 as compared with that without mirror oscillation, and 
the associated ZB frequency should be decreased by a factor 0.7 or 0.5. This is again in agreement with our direct 
numerical experiments. Note that if we have i/o(^^j^) = 0, i.e., exactly on a CDT point [15| where the transition 
between the two dark states is totally suppressed by the driving, then in the case of Fig. 2(a), the ZB amplitude is 
zero, and in the case of Fig. 2(b), the ZB frequency is zero. So in either case the ZB effect disappears on a CDT 
point. Interestingly, this is also a situation where -ffdrivcn takes a spin-helix form [l3|, which now possesses an SU(2) 
symmetry. Such a dynamical realization of a spin-helix Hamiltonian constitutes an intriguing consequence of the 
coupling between the mechanical mirror oscillation and the effective spin-orbit interaction. 

Results in Fig. 2 also indicate that in general cold-atom ZB in a two-dimensional geometry suffers from quick 
damping [l^, thus setting limitations to its potential application. This ZB damping can have two different, but 
related, interpretations. First, Eq. involves an integral over a continuous distribution of the ZB frequency i^(p). 
As time evolves, oscillations associated with different p comprising the initial wavepacket will necessarily dephase. 
This makes it clear that the lifetime of the ZB is inversely proportional to the momentum spread in the initial 
wavepacket The second picture is more enlightening. Note that the oscillating term in Eq. ([6]) is due to the 
quantum interference between the IV"^) and branches in Eq. ([5]). Because these two terms have different group 
velocities Vpi!^±(p, p) = p/m± ?iKVp|p|/m, the wavepacket on the branch will move away from that on the 
\'4>~) branch (see insets in Fig. 2). As a consequence the overlap of the two wavepackets decreases with time, their 
quantum interference decays, and hence ZB damps. The CDT condition j7o(^£^) = may eliminate this group 
velocity difference, but as mentioned above, it also kills ZB in the beginning. 

To dynamically suppress the ZB damping of a two-dimensional cold-atom wavepacket, we next consider a different 
regime of tOd- Suppose the initial wavepacket has an average momentum p^ (p^) along the x (z) direction, with 
pJJ 3> We let LUd satisfy the following "resonance" condition: « mtOd/^K, with ujd 3> Kp^^/m, nvd- In this regime, 
ujd matches the precession frequency of the effective spin around p^. However, in a rotating frame, where the spinor 
wavefunction (a', b')"^ is related to the wavefunction (a, b)"^ in a non-rotating frame by 



1 \ a-b f 1 , ^ 

1 e ^ +^ _i e ^ , (9) 
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FIG. 3: (Color Online) (a) Suppression of ZB damping by mirror oscillation in an on-resonance regime. The initial momentum 
spread and Ud are the same as in Fig. 2. Panels (b) and (c) demonstrate that the wavepacket does not split despite the ZB. 



iOd still dominates over all other frequencies. One can then apply again the above-introduced high-frequency approx- 
imation and obtain the following effective Hamiltonian in the rotating frame (for Vdn/ujd ^ 1) 



TTCS ^ P_ 



Tin 



mvd 



mujd 
2k 



(10) 



on-resonance 

1/2 



Hamiltonian are 



given by £;5?"°(p) = /2m ± 

P^/2to ± hKVd/2, where we have used the resonance condition 



Eigenvalues of this effective 

{Tin/m) {mvd/2f + {jp^ ~ mLJd/2ny 

p° « mujdl2K and the assumption that the Px of a wavepacket is strongly peaked at Same as in our early 
consideration, the beating frequency between the two eigenvalue branches, i.e., w''°''°(p) = E'^^°{p)/'h — E'^_^^°{p)/h, 
gives rise to the new ZB frequency a;'''''^°(p) « nvd- This indicates that now the ZB frequency is totally determined 
by the velocity amplitude Vd of the mirror oscillation, thus entirely converting a mechanical property of the mirror to 
that of a fundamental quantum coherence phenomenon. Further, it is easy to find that (r) here can still be given by 
Eq. (|n]), but with p = {px — mvd/2K)ex + {mvd/2)ez ~ {mvd/2)ez- Further evaluating Vp0(p) in Eq. ([B]), one finds 
that the ZB here should be along the x direction. 



Significantly, the group velocities for both the eigenvalue branches i?5?''°(p) now become the same, i.e 



(P) 



Hence, an initial wavepacket undergoing ZB will not split into two parts. According to our early explanation of the 
ZB damping seen in Fig. 2, such suppression of wavepacket splitting should suppress the damping of ZB. Numerical 
results in Fig. 3 directly using Hdmcn confirm our predictions based on H^l^^^. In particular, the initial state used in 
Fig. 3 has the same momentum spread as in Fig. 2, and the mirror oscillation is now under the resonance condition 
Px « muJd/2K. It is seen from Fig. 3 that the wavepacket does not split despite the ZB. For a similar ZB frequency 
as in Fig. 2(a), the ZB damping can be hardly seen in Fig. 3, even after doubling the time scale. It can be estimated 
that via the mirror oscillation the ZB lifetime here is increased by more than one order of magnitude. 



IV. CONCLUSION 



To conclude, by considering mirror oscillation in a laser-atom system, we have proposed to explore time-dependent 
Dirac-like equations via driving an effective spin-orbit interaction. Using cold-atom ZB as a case study, we have shown 
how mirror oscillation can be used to control the amplitude, the frequency, and the lifetime of ZB in a two-dimensional 
geometry. For m ~ 10^^^ Kg, n ^ 10^ m~^, |p| ^ 1 — 10 we find ^ 10** Hz, falling in the range of the mirror 
oscillation frequency in current optomechanical systems [IQ] . This also suggests that the peak velocity of the mirror 
is of the order of 10~"^ — 10~^ ms~^. In addition to opening up a new means of matter-wave manipulation, our 
theoretical work should greatly motivate experimental efforts in realizing and exploring cold-atom Zitterbewegung. 
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